SPECTRAL PROPERTIES OF 4-DIMENSIONAL 
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Abstract. We study the spectral properties of a large class of compact flat Riemannian man- 
ifolds of dimension 4, namely, those whose corresponding Bieberbach groups have the canonical 
lattice as translation lattice. By using the explicit expression of the heat trace of the Laplacian 
acting on p-forms, we determine all p-isospectral and L-isospectral pairs and we show that in 
this class of manifolds, isospectrality on functions and isospectrality on p-forms for all values of 
p are equivalent to each other. The list shows for any p, 1 < p < 3, many p-isospectral pairs 
that are not isospectral on functions and have different lengths of closed geodesies. We also 
determine all length isospectral pairs (i.e. with the same length multiplicities), showing that 
there are two weak length isospectral pairs that are not length isospectral, and many pairs, 
p-isospectral for all p and not length isospectral. 



Introduction 

If M is a compact Riemannian n-manifold, let Spec p (M) denote the p-spectrum of M, that 
is, the collection of eigenvalues, counted with multiplicities, of the Laplacian acting on smooth 
p-forms. If Spec p (M) = Spec p (M'), then M and M' are said to be p-isospectral. In the function 
case, i.e. if p = 0, one says that M and M' are isospectral. The set of all lengths of closed geodesies 
of M (resp. lengths with multiplicities) is called the weak length spectrum or L-spectrum (resp. 
length spectrum or [L]- spectrum) of M. 

In [6, 7] multiplicity formulas for the eigenvalues of the Laplacian acting on natural vector 
bundles over compact flat manifolds were given, together with isospectrality criteria, that allow 
to give a variety of examples and counterexamples. In the multiplicity formulas for eigenvalues 
for the Laplacian acting on p-forms, certain traces appear. For flat manifolds of diagonal type 
(see Section 1) these are given by integral values of Krawtchouk polynomials, hence one can use 
their integral zeros to construct examples of p-isospectral manifolds that are not 0-isospectral. 
In [8] different types of length spectra are studied and p-isospectrality is compared with length 
isospectrality. Some of the examples in [6, 7, 8] occur already in dimension 4 and, given that 
there is a full classification of the corresponding Bieberbach groups in this case ([1]), it seems 
natural to study the spectral properties of flat manifolds for n = 4. This is the main goal of the 
present paper. 

It is known that if two flat manifolds are isospectral, the corresponding covering tori must 
be isospectral to each other ([9]). It is thus natural, when comparing spectra, to consider 
flat manifolds both having the same covering torus. In this paper we will restrict ourselves to 
manifolds covered by the standard torus, that is, the corresponding Bieberbach groups have the 
canonical (cubic) lattice of translations. We will be left with 54 out of the 74 diffeomorphism 
classes of four dimensional manifolds and we shall see that some of the diffeomorphism classes 
allow several different isometry classes. If we include representatives for all these isometry 
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classes, then we end with a total of 79 different isometry classes of Bieberbach 4-manifolds, 
covered by the standard 4-dimensional torus. 

An outline of the paper is as follows. In Section 1 we briefly review the main facts on 
Bieberbach groups and some results from [6, 7, 8]. We also give a table with the values of the 
Krawtchouk polynomials K^{j) for < p,j < 4 (see (1.3)). 

In Section 2 we give a full set of representatives for the isometry classes in each diffeomorphism 
class, and compute the ingredients for the p-heat traces and the lengths of closed geodesies for 
all Bieberbach groups in our class. At the end of the section, we use formula (1.1) to give explicit 
computations of multiplicities of eigenvalues in the particular case of two flat manifolds with 
holonomy group D4. These turn out to be p- isospectral if and only if p = 1,3. Furthermore they 
do not have the same lengths of closed geodesies. This is the first example with these properties 
for flat manifolds having nonabelian holonomy groups. 

In Section 3 we give expressions for the p-heat traces (s) in the context of this paper and 
then rewrite Zp(s) in terms of polynomial expressions that are simpler than the zeta functions, 
and still encode the same spectral information. 

In the last section we give all pairs of Sunada-isospectral, p-isospectral and L-isospectral flat 
manifolds of dimension 4 in our context, by comparison of the zeta functions Zp(s). We will 
show that all O-isospectral pairs are actually L-isospectral, p-isospectral for all p and have the 
same holonomy representation. Moreover, the non-diffeomorphic pairs occur only in the class 
of Bieberbach manifolds of diagonal type. On the other hand, our list will reveal a large num- 
ber of p-isospectral pairs for either p = 1 and p = 3, or for p = 2, that are not isospectral 
on functions and have different lengths of closed geodesies. Also, we shall see that there exist 
several L-isospectral pairs that are not isospectral. Finally, we take into account length mul- 
tiplicities, finding all [L]-isospectral pairs, i.e. having the same lengths of closed geodesies and 
the same multiplicities for each length. The list obtained shows, already in dimension 4, that 
non diffeomorphic Bieberbach manifolds that are p-isospectral for all < p < 4, are often not 
[L] -isospectral (4 pairs out of 9). 



1. Preliminaries 

We shall first recall some standard facts on flat Riemannian manifolds (see [3]). A discrete, 
cocompact subgroup T of the isometry group of R n , J(IR n ), is called a crystallographic group. If 
furthermore, T is torsion-free, then T is said to be a Bieberbach group. Such T acts properly dis- 
continuously on M n , thus Mp = T\IR n is a compact flat Riemannian manifold with fundamental 
group T. Any such manifold arises in this way. Any element 7 £ I(M. n ) decomposes uniquely as 
7 = BLf,, with B 6 O(n), b 6 M. n and denotes translation by b. The translations in T form 
a normal maximal abelian subgroup of finite index, identified with A, a lattice in M n which is 
5-stable for each BL^ £ T. The quotient F := A\T is called the holonomy group of T and gives 
the linear holonomy group of the Riemannian manifold Mp. The action of F on A defines an 
integral representation of F, usually called the holonomy representation. 

If BLf, is in T, denote rig := dimker(i? — Id). The torsion- free condition on T implies that 
for any BL^ G T we have ub > and furthermore, b + := PB(b) / 0, where ps denotes the 
orthogonal projection onto ker(f? — Id). 

We now recall from [6, 7] some facts on the spectrum of Laplace operators on vector bundles 
over flat manifolds. If r is an irreducible representation of K = 0(n) and G = 7(R n ) we form the 
vector bundle E T over G/K ~ M n associated to r and consider the corresponding bundle T\E T 
over ryiR™ = Mp. Let — A T be the connection Laplacian on this bundle. For any nonnegative 
real number /x, let A* = {A 6 A* : ||A|| 2 = /i}, where A* denotes the dual lattice of A. In [7], 
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Theorem 2.1, it is shown that the multiplicity of the eigenvalue 47r 2 /i of — A r is given by 

(1.1) ^(r) = ^ trr(B)e M)7 

7=BL b £A\r 

where e^ )7 = X^ueA* -Bv=v e~ 2mv ' b . In the case when r = t p , the p-exterior representation of 
O(n), we shall write tr p (S) and dp tfl (T) in place of trr p (5) and d Tpt ^(T) respectively. If p = 
we have A ro = A, the standard Laplacian on functions. 

For a special class of flat manifolds the terms in formula (1.1) can be made more explicit. 
A Bieberbach group T is said to be of diagonal type if there exists an orthonormal Z-basis 
{ei, . . . , e n } of the lattice A such that for any element BL^ G T, Bei = ±e{ for 1 < i < n (see 
[7], Definition 1.3). Similarly, Mr is said to be of diagonal type, if V is so. We note that it may 
be assumed that the lattice A of T is the canonical lattice. 

These manifolds have, in particular, holonomy group F ~ Z£, for some 1 < r < n — 1. 
After conjugation by a translation we may assume furthermore that 6 6 ^A, for any BL^ G T 
(see [7], Lemma 1.4). In this case, the terms e„ 7 in the multiplicity formula (1.1) become 
sums of l's and — l's. Moreover, the traces tr p (B) are given by integral values of the so called 
Krawtchouk polynomials K™{x) (see [6], Remark 3.6, and also [7]; see [4] for more information 
on Krawtchouk polynomials). Namely, we have: 

p 

(1.2) tr p (B) = K;{n - n B ), where K%{x) := ^(-1)* 

t=o 

For further use, we now give a table with the values of the Krawtchouk polynomials K™(j) for 
0<p,j<4. 



(1.3) 



V 





1 


2 


3 


4 




1 


4 


6 


4 


1 


K 4 p (l) 


1 


2 





-2 


-1 


K 4 P {2) 


1 





-2 





1 


K 4 P (3) 


1 


-2 





2 


-1 




1 


-4 


6 


-4 


1 



The next lemma gives some auxiliary facts on the Krawtchouk polynomials. 

Lemma 1.1 (See [4], p. 76). If 1 < k,j < n we have 

(i) K T k l (j) = (-l) j K£_ k (j) = (-l) k K%(n - j). Hence if n even, then K2(j) = for j odd 
and #£(§) = 0, fork odd. 

(ii) (p^(j') = (l)K?(k). This implies that K%(j) = if and only if K?(k) = 0. 

If T is of diagonal type and BLb £ T, we have ub = dim(M n ) B = |{1 < i < n : Bei = ej}|. 
We set 

(1.4) n B{^) := \{l < i < n : Bei = &i and b ■ ei = \ mod Z}|. 
If < t < d < n, the Sunada numbers for T are defined by 

(1.5) c dit (r) := |{BL fe G r : ng = d and n B (I) = t} |. 
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We note that by Lemma 1.1 in [8], one has that ub{\) > for any 7 = BL^ G T, 7 A. 

Remark 1.2. In [7] it is shown that the equality of the Sunada numbers c^(r) = Cd,t(T') for 
every d, t, is equivalent to have that Mp and Mp verify the conditions in Sunada's theorem. In 
this case one says that Mr and Mp are Sunada isospectral (see [7], Definition 3.2, Theorem 3.3, 
and the discussion following it). In particular this implies that Mp and Mp are p-isospectral 
for all < p < n. 

For T a Bieberbach group, r a finite dimensional representation of O(n) and Re(s) > 0, we 
consider the heat trace zeta function 

(1.6) Z?(s) := £ e- As = £ d T , M (r) e" 4 ^. 

AeSpec T (M) /i>0 

This series is uniformly convergent for i?e(s) > e, for any e > 0. In the case when r = t p , for 
some < p < n, we write Zp(s) = Z^ p (s). 

The lengths of closed geodesies in Mr are the numbers + A+|| with BL^ £ T running 
through a full set of representatives of F = A\T and A £ A. Recall that b + = pbQ>), where ps 
denotes the orthogonal projection onto ker(B — Id). 

We will make use of the next result from [8] . 

Theorem 1.3. (i) In the notation above, we have 



|A 1 +6. 



ryY, \ 1 trr(B) us 

BL b eF TO H A ^ A + Gp B (A) 

Spec (Mr) determines the lengths of closed geodesies of Mr and the numbers ng; Spec T (Mr) 
(in particular Spec p (Mr) , for any fixed p > 0) determines the spectrum of the torus Ta = A\R n 
and the cardinality of F. 

(ii) If r is a Bieberbach group of diagonal type with F ~ Z r 2 then we have: 



1 n a 

z p(s) = VrzZ K v( n ~ d ) (H^EQ, t (r)^ t (i) 

where, for Re(s) > 0, 



2 r 

d=l t=0 



(m 1 ,...,m d )eZ d 

Furthermore, M-p and Mr' are p-isospectral if and only if 

k;(u - d) Cd>t (T) = k;{u - d) c dtt {V) 

for each 1 < t < d < n. In particular, if Q jt (r) = c d t (T') for every d,t, then Mr and Mp' 
are p-isospectral for all p. If, for p fixed, K r p l {x) has no integral roots and Mr and Mr' are p- 
isospectral then they are Sunada isospectral, hence q-isospectral for all < q < n. In particular, 
for groups of diagonal type, isospectral implies Sunada isospectral. 

Remark 1.4. In [8] it is shown by asymptotic methods that the formula in (i) of the theorem 
implies that two isospectral flat manifolds must be both orientable or both nonorientable. This 
is not true for p-isospectral pairs (see [6] ) . 
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2. BlEBERBACH GROUPS WITH CANONICAL TRANSLATION LATTICE 

In [1] a complete list of isomorphism classes of the crystallographic groups of dimension 4 
is given. This list contains 4783 groups out of which 74 are Bieberbach groups. Among these 
Bieberbach groups there are 54 which allow the canonical lattice as lattice of translations and 
20 that do not. 

In this section we shall briefly explain how to obtain the Bieberbach groups from the data 
given in Tables 1C and 2C in [1] and how to present them in the notation given in Section 
1. These will be given in tables containing the matrices Bi and the corresponding translation 
vectors bi and some additional information. For convenience, in the last table we give the Sunada 
numbers for those groups T having diagonal holonomy representation. 

For simplicity, we shall denote by 1, . . . ,74, the 74 Bieberbach groups in Table 1C, p. 408 
of [1], giving a full set of representatives for the isomorphism classes. Groups having the 
same holonomy representation (i.e. the same Z-class) are put together into families J-\, ■ ■ . ,J~8- 
In the notation in [1], the Z-class is indicated by the first 3 numbers, for instance, in the 
group 2=2/1/1/2 below, by 2/1/1. The full list is the following: 1=1/1/1/1; 2=2/1/1/2; 
3=2/1/2/2; 4=2/2/1/2; 5=3/1/1/2; 6=3/1/2/2; Fi={7, 8, 9, 10, ll}={4/l/l/6, 7, 10, 11, 13}; 
12 = 4/1/2/4; T 2 = {13,14,15} = {4/1/3/4,11,12}; 16 = 4/1/4/5; 17 = 4/1/6/4; F z = 
{18,19,20,21} = {4/2/1/8,11,12,16}; 22=4/2/3/4; 23 = 4/3/1/6; T± = {24,25,26,27} = 
{5/1/2/7,8,9,10}; 28 = 5/1/3/6; 29 = 5/1/4/6; 30 = 5/1/6/6; 31 = 5/1/7/4; 32 = 5/1/10/4; 
T*> = {33,34,35,36,37,38,39,40,41,42}={6/1/1/41,45,49,63, 64,66,81,82,83,92}; T % = 
{43,44}={6/2/l/27,50}; 45=7/2/1/2; 46=7/2/2/2; 47=8/1/1/2; 48=8/1/2/2; 49=9/1/1/2; 
50=12/1/2/2; 51=12/1/3/2; 52=12/1/4/2; 53=12/1/6/2; 54 = 12/3/4/6; 55=12/3/10/5; 
56=12/4/3/11; F 7 ={57, 58}={13/l/l/8, 11}; 59=13/1/3/8; ^ 8 ={60, 61, 62}={13/4/l/14, 20, 
23}; 63=13/4/4/11; 64=14/1/1/2; 65=14/1/3/2; 66=14/2/3/2; 67=14/3/1/4; 68=14/3/5/4; 
69=14/3/6/4; 70=15/1/1/10; 71=15/4/1/10; 72=24/1/2/4; 73=24/1/4/4; 74=25/1/1/10. 

In Table 1C in [1] we find, for each group V, matrices Ai,...,Ak and translation vectors 
ai, . . . , Ofc together with a (Bravais) type of lattice. The matrices are not orthogonal in general, 
however they are orthogonal with respect to some suitable inner product. In Table 2C in [1], 
there is a list of real symmetric matrices G = (gij) (Gram matrices) corresponding to the types 
of lattices given in Table 1C. The entries of these matrices G are parameters representing the 
scalar products of the Z-basis vectors of a lattice of translations A of any group having this 
lattice type. 

More precisely, Table 1C gives a crystallographic group T\ = {A±L ai , . . . , A k L ak , L^) where 
Ai = 7Lv\ © • • • © 7Lv± is a lattice in R 4 with AiA\ = A\ for 1 < i < 4. As mentioned before, 
Table 2C gives a symmetric matrix G = (gij), where g^ = Wi ■ Wj and A2 = 7Lw\ © • • • © TLw^ 
that gives rise to a group T2 = XTiX^ 1 for some unimodular matrix X. Such X G SL^(Z) 
and its inverse X~ l are given in Table 2C, together with the matrix G. Now, the matrices 
B 1 := XAiX- 1 , ...,B k := XA k X~ x obtained in this way become orthogonal with respect to 
the basis {wi, . . . , W4} of A2. We will take A2 = Z 4 . We list the matrices X / Id at the end of 
the section. 

We remark that there is often a difference between the vectors aj in Table 1C and those used in 
this paper. This is so because, as explained in Section 1, any element 7 £ T decomposes uniquely 
as 7 = BLf,, with B G O(n) and b G W\ and so 7 acts on x € M. n by j(x) = Bx + Bb. However, 
in [1] the elements of V are written in the form 7 = (A, a) with action j(x) = Ax + a, so (A, a) 
corresponds to AL^-i a in our notation. We still have to conjugate (A, a) by the corresponding 
matrix X. In this way, from the original data (A, a) we obtain BLf, = XAX~ l L X A- 1 a- 

In this paper we are interested in groups having the canonical lattice as lattice of translations. 
Therefore, we need to look for the types of lattices having Gram matrices G whose parameters 
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admit the identity matrix G = Id. For example, for the lattices of type X/IV (see [1], p. 275) 
we have 



and we see that G ^ Id for any possible choice of a, b, c. Thus, a Bieberbach group having type 
X/IV does not admit the canonical lattice as translation lattice, and hence the groups 52, 59, 
63 are not in the class considered. 

Taking into account the full list of Gram matrices in Table 2C in [1] one can check that the 
4-dimensional Bieberbach groups admitting the canonical lattice as lattice of translations are: 



1, 2, 3, 4, 5, 6, = {7,8,9,10,11}, 12, T 2 = {13,14,15}, F 3 = {18,19,20,21}, 22, 23, 
T± = {24,25,26,27}, 28, 29, F 5 = {33,34,35,36,37,38,39,40,41,42}, F 6 = {43,44}, 45, 
47, 50, 51, 54, 56, JF 7 = {57,58}, F 8 = {60,61,62}, 64, 67, 72 and 74. 



We shall present the above 54 groups in tables, together with 25 additional groups to give 
a complete set of representatives of the isometry classes of groups in our class. Each one of 
the groups in question lies in one table and groups having the same holonomy representation 
are put together in the same table. Also, for convenience, all manifolds having holonomy group 
isomorphic to %2 are put together in the same table (Table 1). Finally, the groups numbered 72 
and 74 are not included in the tables, for simplicity. They can not come up in any p-isospectral 
pair because they are the only groups with holonomy groups of order 12 and 24 respectively. The 
tables show the non trivial elements in the holonomy group, together with the corresponding 
translation vectors and some additional information that will allow to obtain explicit expressions 
for the zeta functions, and to determine the L-spectra of the associated manifolds. 

The tables are organized in the following manner. Suppose that V = (BiL^, . . . , B^L^, A). 
In the first row of the table, labeled £>, we give the non trivial matrix elements of the holonomy 
group F. These matrices B £ 0(4, M.) n GL^Z) are written in the form d(Ci, C2), d{C\, C2, C3) 
or d(Ci, C2, C3, C4) where, for instance, d(C\, C2) indicates a matrix with Ci, C2 lying diagonally 
in B and similarly in the other cases. Here Cj, 1 < i < 4, will be one of the following matrices: 
±1, ±1, ±1, ±J, ± J, ±T, ±T t or K with 



The second group of entries, corresponds to the translation vectors b. Here we include one 
row for each group in the family, denoting with primes the different isometry classes in the same 
diffeomorphism class. 

In the table we also indicate a set of generators over Z of the space fixed by B, i.e. A B = 
{v £ A : Bv = v}, the numbers ub, vol(A B ) and the p-traces tr p (B). In general, these traces 
are given by integer values of Krawtchouk polynomials, Kp(j), but there are some exceptions. 
In these cases all the values of the p-traces, for < p < 4, are listed. The components of the 
vectors b that are not in A B are put in between parentheses in the tables. In the next section we 
will use these tables to compute, for each group T, the corresponding p-heat trace zeta function 
from which one can read off spectral information on Mp. 

Given a group k, there may be several different isometry classes of Bieberbach groups in the 
isomorphism class of k, to be denoted by k', k" and so on (see for instance, Table 2.1). One 
shows that all such classes can be obtained by conjugation of the matrices in the point group F 
of k, by some C £ GL4(Z) n Z(F) \ 0(4), where Z(F) denotes the centralizer of F. 
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The following is a list of representatives for the isometry classes of Bieberbach groups having 
canonical lattice of translations. The corresponding classes in the notation in [1] are given in 
between parentheses. 



Table 2.1: 2 (2/1/1), 3 (2/1/2), 4 (2/2/1), 5 (3/1/1), 6 (3/1/2) (F ~ Z 2 ). 



# 


2 


3 


4 


5 


6 


R 

1 J 


rMT 1 —1) 


r1(T H 


r/f— T — 1 1"! 


rl(-T T) 




b 


e3 
2 


ei 
2 


e4 
2 


e4 
2 


ei 
2 


V 


e2+e3 
2 


ei+e2 
2 




e3+e4 
2 


ei+e 2 +e3 
2 


b" 


ei+e 2 +e 3 
2 


ei+e 3 +e 4 
2 








b'" 




ei+e2+e3+e4 
2 








A B 


ei,e 2 ,e 3 


ei,e 2 ,e 3 + e 4 


e 4 


e3,e 4 


ei,e 2 + e 3 


riB 


3 


3 


1 


2 


2 


vol{A B ) 


1 


V2 


1 


1 


V2 


tr p {B) 




^ P 4 (l) 






«2(2) 


Pi, 02 


3,3 


3,3 


1,3 


2,2 


2,2 


Orientable 


no 


no 


no 


yes 


yes 


Diagonal-type 


yes 


no 


yes 


yes 


no 



Table 2.2: 47(8/1/1) (F ~ Z 3 , non-diagonal, orientable, ft = 2, ft = 1) 





7i 


7? 




d(l,T) 


d(l,T*) 


647' 


ei 
3 

ei+e2+es+e 4 
3 


9 ei 

Z 3 
ei+e2+es+e4 
Z 3 


A B 
n B 
vol(A B ) 


ei, e 2 + e 3 + e 4 

2 

/q 
V ^ 


ei , e 2 + e 3 + e 4 
2 

v/3 


tr p {B) 


110 11 


110 11 
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Table 2.3: Tx = {7, 8, 9, 10, 11} (4/1/1) (F ~ diagonal, non-orientable, Pi = 2, (3 2 = 1) 





7i 


72 


7i 72 


B 


d(/,l,-l) 




d(I,-I) 


bv 

u 7 


2 


e2 
2 


e2+(e 3 ) 
2 




ea 
2 


ei+e 2 
2 


ei+e 2 +(e 3 ) 
2 


/)o 

°8 


e3 
2 


e2+e4 
2 


e2 + (e3+e4) 
2 


bg/ 


e3 
2 


ei+e2+e4 
2 


ei+e 2 +(e3+e 4 ) 
2 


b 9 


£2 
2 


£l 

2 


ei+e 2 
2 


b 9 , 


e2 
2 


ei+e 2 
2 


ei 
2 


ho 


e2 
2 


ei+e4 
2 


ei+e2 + (e 4 ) 
2 


bw 


e2 
2 


ei+e 2 +e 4 
2 


ei+(e 4 ) 
2 


t. 

010" 


ei+e2 
2 


ei+e4 
2 


e2+(e4) 
2 


011 


e2+e3 
2 


ei+e4 
2 


ei+e2+(e 3 +e 4 ) 
2 


hv 


e2+e3 
2 


ei+e2+e4 
2 


ei + (e 3 +e 4 ) 
2 


A B 


ei,e 2 ,e 3 


ei,e 2 ,e 4 


ei,e 2 


n B 


3 


3 


2 


vol(A B ) 


1 


1 


1 


tr p {B) 


K 4 p (l) 


^ P 4 (l) 


K 4 P (2) 



Table 2.4: 12 (4/1/2) (F ~ Z|, non-diagonal, non-orientable, 0i = 2,/3 2 = 1) 





7i 


72 


7i 72 




d(J, J) 


d(I,-J) 


d(J,-J) 


&12 
&12' 


e2 
2 

ei+e 2 
2 


ei 
2 
ei 
2 


ei+e2 
2 

e2 
2 


A B 
wZ(A B ) 


ei,e 2 ,e 3 + e 4 
3 

V2 


ei,e 2 ,e 3 - e 4 
3 

V2 


ei,e 2 
2 

1 


*r p (B) 


^(1) 


K 4 p (l) 
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Table 2.5: T 2 = {13, 14, 15} (4/1/3) (F ~ 1% non-diagonal, non-orientable, fa = l,fa = 0) 





7i 


72 


7i 72 


B 


d(U,-l) 


J, 1) 


d(l, J,-l) 


bis 


2 


(34 

2 


ei+(e4) 
2 


b-i q/ 


ei+e 2 +e 3 
2 


e 4 
2 


ei+e2+e3+(e 4 ) 
2 


7 

014 


e2+e3 
2 


ei 
2 


ei+e2+e3 
2 


bi4' 


e2+e3 
2 


ei+e2+e 3 
2 


ei 
2 


bis 


e2+e3 


ei+e4 


ei+e2+e3+(e4) 

o 
z 


bis* 


e2+e3 


ei+e2+e3+e4 


ei+(e4) 




2 


2 


2 


A B 


ei,e 2 ,e 3 


ei,e 2 + e 3 ,e 4 


ei,e 2 + e 3 


n B 


3 


3 


2 


vol(A B ) 


1 


V2 




tr p (B) 




K${1) 


«2(2) 



Table 2.6: JF 3 = {18,19,20,21} (4/2/1) (F ~ Z|, diagonal, non-orientable, fa = 1, fa = 1) 





7i 


72 


7i72 




d(/,-l,l) 


d(-/,-l,l) 


d(-/,J) 


&18 

&19 
&19' 
&20 
^20' 
&21 
&21' 


e4 
2 

e 2 
2 

ei+e2 
2 

e 2 
2 

ei+e 2 


(e 3 )+e 4 
2 

e4 
2 

e4 
2 

(e3)+e4 
2 

(e3)+e4 


e3 
2 

(e 2 )+e 4 
2 

(ei+e2)+e4 
2 

(e2)+e3+e4 
2 

(ei+e2)+e3+e4 


2 

e2+e 4 
2 

ei+e2+e4 
2 


2 

(e 3 )+e 4 
2 

(e 3 )+e 4 
2 


2 

(e2)+e3 
2 

(ei+e 2 )+e 3 
2 


A B 
n B 
vol(A B ) 


ei,e 2 ,e 4 
3 
1 


e 4 
1 
1 


e 3 ,e 4 
2 
1 


tr p (B) 




K^) 
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Table 2.7: 22 (4/2/3) (F ~ Z|, non-diagonal, non-orientable, fa =l,/3 2 = 0) 







T9 
I £ 


1 1 1 z 




d(l, J,l) 


d(-l,-/,,l) 


d(-l,-J,l) 




ei 
2 


e 4 
2 


(ei)+e4 
2 


622' 


e l+ e 2+ e 3 
2 


£4 
2 


2 


A B 


ei,e 2 + e 3 ,e 4 


e 4 


e2 - e 3 ,e 4 




3 


1 


2 


vol(A B ) 


V2 


1 




tr p (B) 




^(3) 


^(2) 



Table 2.8: 23 (4/3/1) (F ~ Z|, diagonal, non-orientable, fa =0,fa = 1) 





7i 


72 


7i72 




d(-/,l,-l) 


d(-/,-l,l) 


d(J,-J) 


623 


e3 
2 


(e2)+e4 
2 


e2+(e3+e4) 
2 


623' 


e3 
2 


(ei+e2)+e4 
2 


ei+e2 + (e3+e4) 
2 


A B 


e3 


e 4 


ei,e 2 




1 


1 


2 


wZ(A B ) 


1 


1 


1 


*r p (B) 


^ P 4 (3) 


^(3) 


^ P 4 (2) 



TaWe 2.9: JF 4 = {24, 25, 26, 27} (5/1/2) (F ~ Z|, diagonal, orientable, fa = 1,02 = 0) 





7i 


72 


7i72 


B 


d(-l, -1,1,1) 


d(l, -1,-1,1) 


d(-l, 1,-1,1) 


624 


e 4 
2 


(e2)+e4 
2 


e2 
2 


625 


e 4 
2 


ei+(e2) 
2 


(ei)+e2+e4 
2 


626 


es. 
2 


ei+(e 2 ) 
2 


(ei)+e 2 +(e 3 ) 
2 


627 


ea 
2 


ei+(e 2 )+e 4 
2 


(ei)+e2 + (e 3 )+e 4 
2 


A B 


e3,e 4 


e 4 , e 4 


e2,e 4 


n B 


2 


2 


2 


vol(A B ) 


1 


1 


1 


tr p (B) 


^(2) 


*J(2) 


*J(2) 
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Table 2.10: 28 (5/1/3) (F ~ Z|, non-diagonal, orientable, Pi = 1,(3 2 = 0) 



li 





7i 


72 


7i72 


B 

028 


d(l,-J,-l) 

ei 
2 


d(l,J,-l) 

ei+(e 4 ) 
2 


d(l,-/,l) 

e 4 
2 


A B 

n B 
vol(A B ) 


ei, e 2 - e 3 
2 

V2 


ei, e 2 + e 3 
2 

V2 


ei, e 4 
2 

1 


tr p (B) 


K 4 P (2) 


K(2) 


K*(2) 



Table 2.11: 29 (5/1/4) (F ~ Z|, non-diagonal, orientable, fa = 1,(3 2 = 0) 





7i 


72 


7i 72 


B 


d(-l,J,l) 


d(l,J,-l) 


d(-l, 1,1,-1) 


&29 


(-e 2 +e 3 )+e 4 
2 


ei 
2 


(ei)+e 2 +e 3 +(e 4 ) 
2 


^29' 


e3+e4 
2 


ei 
2 


(ei)+C3 + (e 4 ) 
2 


A B 


e2 + e 3 ,e 4 


ei,e 2 + e 3 


e2,e 3 




2 


2 


2 


voZ(A B ) 


V2 


V2 


1 


tr p (B) 


K 4 P (2) 


^(2) 


^(2) 



Table 2.12: 45 (7/2/1) (F ~ Z 4 , non-diagonal, orientable, fa = 2,fa = 1) 





7i 


7i 2 


7? 


B 


d(J,-J) 


d(J,-J) 


d(I, J) 


645/ 


e 2 
4 

ei+e2 
4 


e2 
2 

ei+e2 
2 


3e 2 
4 

3ei+3e2 
4 


A B 
n B 
vol{A B ) 


ei,e 2 
2 

1 


ei, e 2 
2 

1 


ei,e 2 
2 
1 


tr p (B) 


1 2 2 2 1 


^ P 4 (2) 


1 2 2 2 1 
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Table 2.13: 50 (12/1/2)(F ~ Z4, non-diagonal, non-orientable, (3\ = l,/? 2 = 0) 





7i 


«,2 

7i 


„,3 

7i 




d(-l,l,-J) 


d(I-I) 


d(-l,l,J) 




e2 
4 


e2 
2 


3e 2 
4 


A B 


e2 


ei,e 2 


e2 


n B 


1 


2 


1 


vol(A B ) 


1 


1 


1 


tr p (B) 


1 -1 


^(2) 


1 -1 



Table 2.14: 51 (12/l/3)(i ? ~ Z4, non-diagonal, non-orientable, Pi = 0,/? 2 = 1) 





7i 


7i 2 


7? 


&51 


d{J,-J) 

e2 
2 


d(I,-I) 

ei+e 2 
2 


d( J, J) 

ei 
2 


A B 
voZ(A B ) 


ei + e 2 
1 


ei, e 2 
2 

1 


ei + e 2 
1 

V2 


*r p (B) 


1 -1 


^(2) 


1 -1 



Table 2.15: 64 (14/1/1)(F ~ Zg, non-diagonal, non-orientable, j3\ = l,/3 2 = 0) 





7i 


7? 


7? 


7? 


7f 


B 


d(l,-T) 


d(l,T*) 


d(l,-l,-l,-l) 


d(i,r) 


d(i,-r*) 


064 


ei 
6 


ei 
3 


ei 
2 


2ei 
3 


5ei 
6 


A B 


ei 


ei, e 2 +e 3 + e 4 


ei 


ei, e 2 +e 3 + e 4 


ei 


n B 


1 


2 


1 


2 


1 


vol{A B ) 


1 




1 


Va- 


1 


tr p (B) 


110-11 


110 11 


^ 4 (3) 


il 1 1 


110-11 



Table 2.16: 67 (14/3/l)(F ~ D 3 , non-diagonal, orientable, fa = 0,(h = 0) 





7i 


il 


72 


7i 72 


7?72 


B 

b67 


d(l,T) 

ei 1 e3+e4 
3 1 2 


d(i,r*) 

2ei I e2+C3 
3 1 2 


d(-i, J,i) 

e4 
2 


4-1,1, J) 
(f) + f 


4-1, K) 

^2ei-j | e2+e3+e4 


A B 
n B 
vol(A B ) 


ei, e 2 +e 3 +e 4 
2 


ei, e 2 +e3+e4 
2 

Va- 


e2 + e3,e 4 
2 

V2 


e2, e 3 + e 4 
2 

V2 


e 3 , e 2 + e 4 
2 

V2 


tr p {B) 


11011 


il 1 1 


^ P 4 (2) 


^ P 4 (2) 


^ 4 (2) 
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Table 2.17: T h = {33, 34, 35, 36, 37, 38, 39, 40, 41, 42} (6/1/1) (F ~ Zf, diagonal, non-orient- 
able, fa = l,fa = 0). 





7i 


72 


73 


7i 72 


7173 


7273 


717273 




d(I,I) 


d(I,I) 


a!(-/,/) 


d(I,I) 


d(-I,I) 


d(-I,I) 


d(-/,J) 


633 


e 4 




(ei)+e 4 


e2+e4 


ei 


(ei+e 2 ) + e 4 


ei+(e 2 ) 


2 


2 


2 


2 


2 


2 


2 


634 


e 4 
2 


E2 
2 


(ei) + E3+e4 
2 


e2+e4 
2 


ei+e 3 
2 


(e!+e 2 +e3) + e4 
2 


ei+(e2+e3> 
2 


635 


£4 


e2+e 4 


(«l) + e3 


£2. 


ei+e3+e 4 


(ei+e2+e 3 ) + e 4 


ei+(e2+e3> 


2 


2 


2 


2 


2 


2 


2 


636 


e.3 
2 


e 2 
2 


e 4 
2 


e 2 + ( e 3) 
2 


e 3 + e4 
2 


( e 2) + e 4 
2 


(e2 + e 3) + e 4 
2 


637 


e 3 


E2 


(ei)+e 4 


e 2 + (e3) 


ei + e 3 +e 4 


(ei+e 2 ) + e 4 


ei+(e2+e3) + e 4 


2 


2 


2 


2 


2 


2 


2 


638 


e3 


e2 


(ei) + e3+e 4 


e 2 + (e 3 ) 


ei+e 4 


( e l+ e 2+ e 3) + e 4 


ei+(e 2 ) + e4 


2 


2 


2 




2 


2 




639 


e3 


ei+e 2 


e 4 


(ei) + e 2 + (e 3 ) 


«3 + e4 


(ei+e 2 ) + e 4 


ei+(e 2 +e 3 ) + e 4 


2 


2 


2 


2 


2 


2 


2 


640 


e3 
2 


ei+e 2 
2 


(ei) + e 4 
2 


(ei) + e 2 + (e 3 ) 
2 


ei + e 3 +e 4 
2 


(e2)+e4 
2 


(e2 + e3) + e 4 
2 


641 


e3 
2 


ei+82 
2 


(ei) + S3+ e 4 
2 


(ei) + e 2 + (e 3 ) 
2 


ei+e 4 
2 


(e 2 +e 3 )+e 4 
2 


(e2) + e4 
2 


642 


e 3 +£4 


e 2 +e 4 


(ei)+e 3 


e2 + (e 3 ) 


ei+e 4 


( ei +e 2 +e 3 ) + e4 


ei+(e 2 ) 
















A B 


e2, e3, e4 


ei, e2, e4 


e3, e4 


e2, e4 


ei , e3, e4 


e 4 


ei, e 4 


riB 


3 


3 


2 


2 


3 


1 


2 


vol(A B ) 


1 


1 


1 


1 


1 


1 


1 


tr-p(B) 






*2(2) 


^ 4 (2) 




^ P 4 (3) 


^ 4 (2) 



Table 2.18: .F 6 = {43,44} (6/2/1) (F ~ Z|, diagonal, non-orientable, fa = 0,/3 2 = 0) 





71 


72 


73 


7172 


7173 


7273 


717273 


B 


d(/,-/) 


d(I,-I) 


d(-I,I) 


d(I,I) 


d(-I,-I) 


d(-I,-I) 


d(-i,i) 


643 
644 


«3 
2 

e2+e 3 
2 


«2 
2 

e 2 +(e4> 
2 


(ei)+e 4 
2 

(ei)+e 4 
2 


e 2 + (e3> 
2 

(«3) + e4 
2 


(ei) + e3+( e 4> 
2 

(ei+e 2 ) + <33+(e4) 
2 


ei+(e 2 +<!4) 
2 

ei+(e 2 ) 
2 


ei+(e2+e 3 ) + e 4 
2 

ei+(e 3 ) 
2 


A s 
n s 
w>2(A s ) 


ei, e2, ez 
3 
1 


e2 

1 
1 


e3, e4 
2 
1 


e2,e4 
2 
1 


e3 
1 

1 


ei 

1 

1 


ei,e4 
2 

1 


tr p (B) 




^(3) 


^(2) 


^ P 4 (2) 


^(3) 


^(3) 


^ P 4 (2) 



Table 2.19: T-j = 57, 58 (13/1/1) (F ~ Z 2 x Z 4 , non-diagonal, non-orientable, fa = l,fa = 1) 





71 


7? 


7? 


72 


7172 


7? 72 


7i72 


£ 


d(/,-J) 


d(I,-I) 


d(I,J) 


d(l I) 


d(I,-j) 


d(I,-I) 


d(I, J) 


657 


«2 
4 


«2 
2 


3e 2 
4 


e 3 + e 4 
2 


B2 1 f es+e 4 \ 

4 + 1 2 ; 


E2 + (e 3 +e4) 
2 


3e 2 1 ( e 3 + e 4 \ 

4 + 1 2 ; 


658 


«2 


«2 


3e 2 


(-ei) + e3 + e 4 




e2 + (-ei+e 3 +e 4 ) 


3e 2 1 /" ei — e 3 — e 4 "\ 


4 


2 


4 


2 


4 + { 2 J 


2 


4 + l 2 J 


A s 


e2 


ei,e 2 


e2 


e2,e3,e 4 


ei,e 2 


e2 


ei,e 2 




1 


2 


1 


3 


2 


1 


2 


«)«(A S ) 


1 


1 


1 


1 


1 


1 


1 


tr p (B) 


1 -1 


^ 4 (2) 


1 -1 


^(1) 


1 2 2 2 1 


^ 4 (3) 


1 2 2 2 1 
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Table 2.20: 54(12/3/4) (F ~ D 4 , non-diagonal, non-orientable, Pi = 0,/3 2 = 0) 





71 


7? 


7i 3 


72 


7172 


7? 72 


7i72 


B 


d(J, j) 


d(I,-I) 


d(J, — J) 


d( J, -1, 1) 


d(I, J) 


<Z(J, 1, -1) 


<i(7, - J) 




ei + (e 4 ) 




S2 + (e3) 


e 4 




ei+e 2 +e 3 


e l+( e 3+ e 4) 


2 


2 


2 


2 


2 


2 


2 


A B 
n B 
vol(A B ) 


ei + e 2 
1 


ei, e 2 
2 

1 


ei + e 2 
1 


ei + e 2 , e 4 
2 


ei, e 2 , e 3 + e 4 
3 

V2 


ei + e 2 ,e 3 
2 


ei, e 2 , e 3 - e 4 
3 


tr p (B) 


1 -1 


^(2) 


1 -1 


^(2) 




^ P 4 (2) 





Table 2.21: 56(12/4/3) (F ~ D4, non-diagonal, non-orientable, fa = 0,02 = 0) 





71 


7? 


7? 


72 


7172 


7? 72 


7? 72 


f>56 


d(J, J) 


d(J,-J) 

S2 + (e3+e 4 ) 


d(I,-J) 

3|2 


d(-J,-J) 

2 


d(-I,I) 


J) 

ei+(e 2 ) + e 3 + e 4 


d(-I,-I) 

^2e!+3e 2 -j | e 3 


2 




2 




A s 

M)2(A S ) 


e 2 
1 
1 


ei,e 2 
2 
1 


e 2 
1 

1 


ei , e3 — e 4 
2 

^2 


e 4 
1 

1 


ei , e 3 + e 4 
2 


e 3 
1 

1 


tr p (B) 


1 -1 


^P 4 (2) 


1 -1 


^ P 4 (2) 


^(3) 


^P 4 (2) 


^(3) 



Table 2.22: F 8 = {60,61,62} (13/4/1) (F ~ D 4 , non-diagonal, orientable, fa = l,fa = 0) 





71 


7? 


7? 


72 


7172 


7? 72 


7?72 


B 


d(J, J) 


d(I,-I) 


d(I,-J) 


d(I, J) 


d(I,-I) 


d(I,-J) 


d(/,J) 


&60 
t>61 

6a 2 


ei 
4 

) 

ei | e 2 + (e 4 ) 
4 "f" 2 


ei 
2 

e l+( e 3+ e 4) 


4 

3e l / e.3 ~i 
4 2 ^ 

3ei | e 2 + (e 3 ) 
4 ' 2 


E2 
2 

«2 
2 

S2 
2 


(f) + f 


(ei)+e 2 
2 

(ei>+«2+e3+e4 


^3exj | e 2 +e 4 


2 

ei+(e3+e4> 
2 


(f) + f 


2 

(eiH<2+e3+e4 
2 




A s 
*)o«(A s ) 


ei, e 2 
2 
1 


ei, e 2 
2 
1 


ei, e 2 
2 
1 


e 2 ,e 3 + e 4 
2 

V^2 


e 2 , e 3 
2 
1 


e 2 , e 3 — e4 
2 


e 2 ,e4 
2 

1 


tr-p(B) 


1 2 2 2 1 


^ P 4 (2) 


1 2 2 2 1 


^(2) 


^ P 4 (2) 


^ P 4 (2) 


^ P 4 (2) 



Example 2.1 (A 1-isospectral pair with holonomy group D4). As an application of the methods 
in this section, we will use formula (1.1) to carry out explicit calculations of multiplicities of 
eigenvalues for the groups 60 and 61, both having holonomy groups isomorphic to D4. As 
a consequence, we shall see that 60, 61 are 1-isospectral and 3-isospectral but they are not 
p-isospectral for p / 1,3. Since the manifolds associated to the groups 60,61 are orientable, 
1-isospectral and 3-isospectral are equivalent in this case. 
According to formula (1.1) we have: 

(2-1) V(r) = iQ|A^ + | Yl tr p (B)e^(T) 

7 =BL 6 GA\r 

where A M = {A G A : ||A|| 2 = fi}. For the elements 7 = 7^, 72, 7172, 7i72, 7i72j we have that 
tr p {B) = K*{2) (see Table 2.22 above), so for p = 1 or p = 3 these traces vanish (see (1.3)). 
Thus, using that tr p {B{) = tr p {B\) = 2 for p = 1, 3, we get 

(2-2) d 1>lt (T) = d 3 ,,(T) = i|A M | + ±e Mi7l (r) + \e^s(T). 
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In order to check 1-isospectrality it suffices to show that e Mi7l (reo) 
01 )■ We have: 

2niv- 

c 

«eA„ 



e^, 7l (r 6 i) and e^^ao) = 



e„ )7 3(r. 



(Tao)= E 



E 



(r 6 i)= £ 



E 



3 — ^zv-ei — mv-ei 



B\v=v 



IHI 2 = 



Since e ~™' e4 = 1 for -u £ Zei©Ze2, we have that e Mi7l (T6o) = e^i^ei)- Similarly, we get that 
e /lj ^(rgo) = e /i, 7 3(r6i)- Observe that we can read off this fact from the rows of the vectors b in 
Table 2.22, since we know that the orthogonal projection of the vectors in between parentheses 
onto the space fixed by the corresponding matrix B, is zero. In this way we have proved that the 
manifolds Mqo,Mq^ obtained from the groups 60,61 are 1-isospectral and hence 3-isospectral. 

We now show that they can not be 0,2 nor 4-isospectral. Take \i = 1. The formula for the 
multiplicities for p = 0, 2, 4 in this particular case is: 

dp,i = g(J)|Ai| + \tr p {Bi) (ei, 71 + e lj7 s) + \K P A {2) (e 1)7 a + ei i7a + e v 

Computing the corresponding ei i7 's we obtain 



■~r,2 + e l,7?72 + e l,7i72 



)• 



7 


7i 


7? 


7? 


72 


7i 72 


7i72 


7i72 


ei, 7 (r 60 ) 


2 





2 


-2 





-2 





ei, 7 ( r 6i) 


2 





2 


-2 


-4 


-2 


-4 



/,From this, and taking into acount that Ai = {±ei; ±e2] ±e3; ±64}, we get 



d p ,i(r, 



61, 



+ h*rp(Bi) 

[ ' + \tr p {B 1 ) 



- K 4 ( 2 ) 
|K p 4 (2). 



Since K^{2) / for p = 0, 2, 4 (see (1.3)), it is clear that dp^{TQo) 7^ <ip, 1(1^1) and this shows 
that M60 and Mq\ can not be p-isospectral for p = 0, 2, 4. For example, for p = 0, 4, we have 
d p ,i{Teo) = 1 and dp^rgi) = 0, while for p = 2 we have ^2,1^60) = 6 and ^2,1(^61) = 4. 



Appendix. Here we list the matrices X ^ Id and C, corresponding to the groups described in 
the tables in this section. 

/1 o\ /o 1 o\ /0100 

AT 3 =X 12 = 0100 , X 6 = 1001 \ X 13 = X 14 = X 15 = X 22 = 1000 

V0011/ Vooio/ V0001 

1\ / 1 

V — V -I 1 v _ -1 1 









( ! \) 


,ci' = 















28 — ^29 — [ -1 I ) ^-47 — 0-10-1 
0010/ VO-100 

c<i — c 1 ' _|ii i r 1 " _ i 1 1 1 i r 1 ' — r 1 ' — n' — r' 

u 2 — L/ 29 — I 1 I ) L/ 2 — I 1 I i°3 - u 9 — ^lO — u ll 

w// _ w _ [ 11 \ w _ / 1 1 \ (~t) _ (~<) _ _ _ _ Ml 

°3 — °47 - I 1 1 I i u 5 ~ I 1 1 I )°6 - u 13 — L/ 14 ~ °15 _ °22 ~ I 1 1 

Z 1 1 

C/ ril rm r~il f-il f-il f-il f~n f~n / 1 

7 — U 8 — °10 — °12 — U 19 — U 20 — U 21 ~ °23 ~ °45 ~ I 1 



16 



R. J. MIATELLO AND R. A. PODESTA 



Finally, we shall give here a list with the Sunada numbers c^t's for those groups having 
diagonal holonomy representation. We only give the values for 01,1,02,1,02,2,03,1,03,2,03,3, since 
ci,o = C2,o = 63,0 = 0, 04,0 = 1 and 04,1 = 04,2 = 04,3 = 04,4 = for any 4-dimensional Bieberbach 
group of diagonal type. 



1 # 


01,1 


t-2,1 


<-2,2 


03,1 


c 3,2 


c 3,3 


2 











1 








1' 


n 


n 


n 


n 


1 

± 


n 




n 
u 


n 
u 


n 
u 


n 
u 


n 
u 


1 
1 


4 


1 

















5 





1 














5 


U 


u 


1 


U 


L) 


I) 


7 


n 


1 


n 


2 


n 


n 


7' 


n 
u 


n 


1 
1 


1 
1 


1 
1 


n 

u 




n 
\j 


1 


n 

yj 


I 




n 





n 
u 


n 

u 


1 
1 


1 
1 


n 


1 
1 


9 








1 


2 








9' 





1 





1 


1 





10 








1 


1 


1 





10' 





1 





1 





1 


10" 





1 








2 





11 








1 





2 





11' 





1 








1 


1 


18 


1 


1 





1 








19 


1 


1 





1 








19' 


1 


1 








1 





20 


1 





1 


1 








20' 


1 





1 





1 





21 


1 


1 








1 





21' 


1 


1 











1 



r # 


Cl,l 


C2,l 


C2,2 


c 3,l 


C3,2 


C3,3 


23 


2 


1 














00/ 
^0 


9 


n 
u 


1 
1 


n 
u 


n 
u 


n 
u 


24 





3 














25 





2 


1 











26 





3 














27 





1 


2 











33 


1 


2 


1 


3 








34 


1 


1 


2 


2 


1 





35 


1 


3 





1 


1 


1 


36 


1 


3 





2 


1 





37 


1 


2 


1 


2 





1 


38 


1 


1 


2 


2 


1 





39 


1 


2 


1 


1 


2 





40 


1 


3 


3 


1 


1 


1 


41 


1 


2 


1 


1 


2 





42 


1 


3 








3 





43 


3 


2 


1 


1 








44 


3 


3 








1 






3. ZETA FUNCTIONS AND p-HEAT TRACE POLYNOMIALS 

In this section we will show that the p-heat trace Zp(s) of an n-dimensional compact flat 
manifold Mr can be viewed as a polynomial of degree n in a finite number of algebraically 
independent functions x(s), y(s), zi(s), . . . , z m (s), for s > 0. We let the dimension n of Mp be 
arbitrary, until further notice. 

We begin by recalling some facts on the L-spectrum of a Bieberbach manifold T. For 7 = BL^, 
denote by ps ■ M n —> (M. n ) B the orthogonal projection and put b + = Pb(6). The L-spectrum 
of Mp is given by SpecL(Mr) = {\\b+ + A+|| : BLi, G T} with BLi, running through a full 
set of representatives of F and A G A (see Section 1). Note that for the -L-spectrum, only 
the components of b not between parentheses in the tables count. In this section, to each 
element B we will associate a monomial in several variables and will express the corresponding 
zeta functions in terms of these monomials, thus getting polynomial expressions for the p-heat 
traces in which the -L-spectrum information is encoded. We will give a complete list of these 
polynomials from which one can recover the L-spectrum of the corresponding manifold. 

By Theorem 1.3, for a compact flat manifold Mr = T\M. n with translation lattice A and 
holonomy group F we have the following expression for the p-heat traces: 

(3-D ^w = A E S (4 ™ r=f £ ^ 

BL b £A\r ^ > A+e Ps (A) 



2 



and if T is of diagonal type, with F ~ ZJj, we can write: 

d 

■ V K n (n.-d\ (4tts)-V 

2 r 



(3-2) Z p r ( S ) = ^if;(n - a!) (4™)"* £ c d , t (r) fc, t (£) 
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where, for Re(s) > 0, 

(mi ,...,m t j)GZ d 

The theta functions 9d,t have simple expressions in terms of 9$ an d @i-> where 
(3.3) 9 (s) := e lfi (s) = ^ e~ sm \ 9 1 (s) := 9 1A (s) = ^ e -^+™) 2 , for Re(s) > 0. 

We define the functions 

x(s) := and y(s) := 



* Aits V 47rs 

Lemma 3.1. For Re(s) > we have 9d t t(s) = 9^~ t (s) 9\{s) . Furthermore, if a^t £ C and 
J2o<t<d a d,tx(s) d ~ t y(s) t = 0, then a d:t = for any < d < t. 

Proof. Decomposing Z d as Z d ~ l © Z* we get 

M«)= E e- S ^=i ( 5 + ^) 2 ^ e- S ^m,- 2 = ^- t(s) ^ i(s) _ 

(mi,...,m t )6Z* (mi,...,m d _ t )eZ d - t 

To prove the second assertion, we note that, as s ^ +oo, 
Now, if u = l/4s then, as s j 0, 



E a,yx(*ril(»)' = 

Kt<d 

E ^ 2 ^)" /2e 



0<t<d 0<<<ci 



tu 
4 



0<<<d 

Now if not all a^t = 0, let t m be minimal with the property that some ad,t m / and let 
be maximal among the d so that a<i,t m 7^ 0. Multiplying the previous expression by u M ' e 4 
and letting u — ► co we have that all terms, except the one corresponding to d — dj^^t — t m , 
tend to zero. Thus, we get that et(2 M ,t m 2*7r d / 2 = 0, a contradiction. □ 

Proposition 3.2. Let T be an n- dimensional Bieberbach group of diagonal type with canonical 
translation lattice A and holonomy group F. Then Z p (s) 6 p^Z[x(s), y(s)] has degree n, no 

independent term, and the coefficient of x{s) k is (resp. 1) for all k = 1, . . . ,n— 1 (resp. k = n). 

_n ||a|| 2 

Proof. By (3.2) and using the notation Zid(s) := (4ns) 2 J2\eA e 4s we can wr ite 

(n\ n ~ l d a (J_) 

z p r (s) = $ z Id (s) + ^ Y1 K > ~ d ) E c ^( r ) 

d=i t=o (4vrs)2 

But 

Z Id (s) = ^s)-2 £ e ^ = -^L = X (sT 

and similarly, 

(3 4) e dAis) _ °o =x ( 8 )d-t t \t 



18 



R. J. MIATELLO AND R. A. PODESTA 



thus 

n— 1 d 

\F\ Z?(s) = Q x(sY + £ K?(n ~d)J2 <*«(T) *(*)*-* vto*. 

d=i t=o 

Since all numbers ( ra ), Kp{n — d) and c^tT) are integers, we are done. □ 

In the case of groups with holonomy representation that is not of diagonal type, one has 
an analogue of (3.2) but the expression is more complicated. Since A is the canonical lattice, 
A = A*. Furthermore, we have that 1 < n B < n and ng = n if and only if B = Id. Thus, 
adding over 7 = BL^ with fixed = d, we may write: 

(3-5) ^( S ) = ^Z H ( S ) + ^g( 4 -)~ f E vSP)^) 



=1 7 = £L 6 e A\r 

n B = d 



l|A + +b + || 2 



where Z 1 {s) := Ea + g Pb (a) e ~ 

Remark 3.3. Let T a Bieberbach group of even (resp. odd) dimension n with diagonal holonomy 
representation and canonical lattice of translations A. If Mr is the associated compact flat 
manifold then Mr is orientable if an only if (x(s), y(s)) has only even (resp. odd) degree 
terms for all < p < n, that is, c<^(r) = for all d = 1, 3, . . . , n — 1 (resp. d = 0, 2, . . . , n) and 
for all t. 

Indeed, let r : V — > O(n) be the projection r(BLb) = B. It is known that Mr is orientable if 
and only if r(j) G SO(n) for all 7 G T, that is, n B G 2Z (resp. n B G 2Z + 1) for all 5 G F since 
1? is a diagonal matrix for each BL^ G T. Now, using the expression (3.5) with canonical A and 
diagonal matrices, and comparing it with the expression (3.2), we see that: 

£ tr p (B)(47rs)-5 Z 7 r ( S ) = £ K p > - d)c^(r) x^^s)' 

7 = BLf, e A\r l<t<d 
n B = d 

are the terms of degree d, and this says that Mr is orientable if and only if Zp (x(s), y(s)) has 
only even (resp. odd) degree terms. Now since 

Coeiixisf-tyisY) = K^{n - d) c d , t (T) 

and Kp(n — d) 7^ in general, it follows that Mr is orientable if and only if c^t(r) = for each 
d = 1, 3, . . . , n — 1 if n is even (resp. d = 2, 4, . . . , n — 1 if n is odd), and for all t. 

Dimension 4- We now look at the special case of dimension 4. From Section 2 we know that the 
information of the L-spectrum is encoded in the terms Z^(s) of Zp (s). From expression (3.2) 
we just have to look at the flat's, or as we have seen, at the monomials x(s) d ~ t y(s)*. In general, 
however, if T is not of diagonal type, the squared lengths of the closed geodesies, ||A+ + 6+|| 2 , 
have summands different from m 2 or (m + ^) 2 . For Re(s) > and r G < r < 1, we 

define the function 

Ms) := Yl e 1 73 ' 



and we set 



(mi,...,m d )eZ d 



Zd,r( S ) : = 



-(c) 



/47TS 

Remark 3.4. Note that (^^(z) = 0o(^) and X 1(2) = 8i(z) so we could just use the 
instead of the to express the zeta functions as polynomials. 
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Proposition 3.5. Let T be a ^-dimensional Bieberbach group with canonical lattice A and 
holonomy group F. Let K = Z(\/2, \/3). Then Z^{s) is given by a polynomial expression 
with coefficients in K in the variables x(s), y(s), z± tr (s), Z2, r (s), Zs, r {s), where r = k/\F\ and 
< k < \F\ - 1. 

Proof. The proposition follows by case by case verification, directly from the tables in Section 
2. We need only check those groups with non-diagonal holonomy representation. We will carry 
out the verification for just one group to ilustrate the method. 

Take the group numbered 29 (see Table 2.11). For the first element B\ we have b\ = — C2 + 
e 3 + e4 and A Bl has the basis {e2 + e 3 , e^}. Since the vector — e2 + e 3 is orthogonal to C2 + e 3 , 
it appears in the table in between parentheses. Thus (&i)+ = \&a- Take A = m\e\ + • • • + 771464 
in A. The orthogonal projection of A onto the subspace generated by e2 + e 3 is given by: 

((m iei + m 2 e 2 + m 3 e 3 + m 4 e 4 ) • ^)) = {m ^ i] (e ^ 3) , 

thus b 1+ + A+ = («^p)(£2_^) + + 1^ In thig way> for 7i = Bl L bl , we have: 

f (Hii^!i2) 2+ ( m 3 + l)2 

4 9 (*) = E e 4s = E e 1 

A + eA s l (mi,m2,m3)eZ 3 

(! ^ )2 (ma+ , )a 

= E e"^^ J] e"" = fc >0 (£) 

(mi,m 2 )eZ 2 m,3eZ 

Doing the same for 72 = B 2 L b2 and 73 = 7172 we get Z^(s) = ^2,0(35) ^1(47) and ^7 3 9 ( s ) = 
#2,2(47) = ^1(47)- Now, using formula (3.5) we get for this case 

^ 29 00 = $ ZM + (f 4 9 (^) + f + ^ 3 9 (*)) 



T( rf 1 ^ 1 
4 *w + 2 4vrs 



(>/2fc,o(£)*i(£) + *?(£)) 



= $ x( S ) 4 + ^ (V2 y(s)z 2 ,o(s) + y(sf) 



4 ^\ J J t 2 

By repeating this argument for each group of non-diagonal type the proposition follows. □ 



We note that these polynomials carry the L-spectrum information only in the variables x(s), 
y(s), Zi tr (s), 1 < i < 3. From Propositions 3.2 and 3.5 the following result is clear: 

Corollary 3.6. M-p and Mp' are L-isospectral if and only if the nonzero monomials appearing 
in Z r (s) and Z r (s) are the same. In particular, ifT, V' are of diagonal type then: Mr and Mr' 
are L-isospectral if and only if for each fixed l<t<d<3we have 

{(M) : c^(r) / 0} = {(M) : c^(rO / 0}. 

We now give the list of p-heat trace polynomials for the groups in the class considered in this 
paper, except for 72, 74, which have holonomy groups of orders 12 and 24 respectively and can 
not be p-isospectral to any other Bieberbach group. 
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Theorem 3.7. In the notation above we have: 
F ~ {Id} : 

Zl(s) = Qx{s)\ 
F ~ Z 2 : 

2Z*(s) = Q x{sf + K±(l) x(s) 2 y(s). 
2Z p 2 '( S ) = ( 4 )x( S ) 4 + K 4 (l) x{s)y{sf. 
2Zf(s) = Qx(sr + K*(l) y(sf. 

2Z*{s) = Q *(*) 4 + f K 4 (l) x(s)y(s)z 2 , (s) = 2Zf(s). 
2Zf{s) = x( S ) 4 + f K 4 (l) y( S ) 2 z 2 , (*) = 2Zf (a). 
2Z p 4 ( S ) = Qx( S ) 4 + ^ 4 (3)y(s). 
2Z p 5 ( S ) = ( 4 )x( S ) 4 + ^ 4 (2) x(a)y( a ). 
2Z p 5 '( S ) = ( 4 )x( S ) 4 + K 4 (2) y( S ) 2 . 
2Z°(s) = Q x(s) 4 + f i^ 4 (2) y( S )z 2>0 ( S ) = 2Z p 6 '(s). 
F~Z 3 : 

3Zf( S ) = ( 4 ) x( S ) 4 + f tr p (B) (z ltl/3 (s)z 3 , (s) + z 1>2/3 (^3,o(s)) = 3Zf(s). 
F^7L\: 

AZl{s) = Q x( S ) 4 + 2K 4 (1) x( S ) 2 y(s) + K 4 (2) x(s)y(s). 

42? (s) = (p x( S ) 4 + K 4 (l) (z(s) 2 y( S ) + x( S )y( S ) 2 ) + i^ 4 (2) y(s) 2 = 4Z p 10 ( S ). 

4Z p 8 ( S ) = Q x(s) 4 + K 4 (l) (x( S ) 2 y(s) + x( S )y( S ) 2 ) + i^ 4 (2) x( S )y( S ) = 4Z p 9 '( S ). 

4Z p 8 '( S ) = ( 4 ) x(,) 4 + K 4 (l) (x(s)"y(s) + 2/(*) 3 ) + ^ P 4 (2) y(*) 2 . 

4Z»(s) = Q x(,) 4 + 2K 4 (1) x(sfy( S ) + K 4 (2) y(sf . 

AZf{s) = Q x( S ) 4 + K*(l) (x( S ) 2 y( S ) + y( S f) + K 4 (2) x( S )y(s). 

4Z p 10 "( S ) = ( 4 ) x(,) 4 + 2K 4 (1) x(s)y(sf + K 4 (2) *( a )y( a ). 

4Z p "(a) = x(,) 4 + 2^ 4 (1) x(,)y( S ) 2 + i^ 4 (2) y( S ) 2 . 

4Z p 11 '( S ) = ( 4 ) x(,) 4 + K 4 (l) (x(s)y(sf + y(,) 3 ) + K 4 (2) x(s)y(s). 

AZ?{s) = Q x( S ) 4 + V2^ 4 (l) x(s)y(s)z 2fi (s) + *J(2) y(s) 2 . 

4Zf(s) = ( 4 ) x( S ) 4 + f K 4 (l) (y( S ) 2 z 2i0 ( S ) + xWy^^oW) + K±{2) x{s)y{s). 

AZ?{ 8 ) = Q x(s) 4 + K 4 (l) (x( S ) 2 y( S ) + ^x(8)y( S )z 2 , (8)) + f K 4 (2) y(s)z 2fi (s). 
AZ™'( S ) = ( 4 ) x( S ) 4 + K 4 (l) (y(,) 3 + f x( S )y(8)z 2 , (8)) + f^ 4 (2) y(,)z 2 , ( S ). 
4Z p 14 ( S ) = ( 4 ) x(,) 4 + K*{1) x( S ) (y( S f + ^ y(s)z 2 , ( S )) + ^K 4 (2) y(s)z 2fi (s) = AZ^'(s). 
AZ^(s) = ( 4 ) x( S ) 4 + K 4 (l) (x(s)y(sf + ^y( s f Z2fi (s)) + f K 4 (2) y( S )z 2 , ( S ) = 4Z p 15 '( S ). 

4Z p 18 ( S ) = ( 4 ) x( S ) 4 + K 4 (l) x( S ) 2 y( S ) + K 4 (2) x(*)y( a ) + K 4 (3) y( S ) = 4Z p 19 ( S ). 
4Zf'( S ) = ( 4 ) x(,) 4 + K 4 (l) x(,)y( S ) 2 + K 4 (2) x( S )y( S ) + K 4 (3) y(s) = 4Z^( S ). 
4Z™(s) = Q x(,) 4 + K 4 (l) x{sfy{s) + K 4 (2) y( S ) 2 + K 4 (3) y{s). 
±Zf\s) = ( 4 ) x(,) 4 + K}(1) x(8)y(s) 2 + K 4 (2) y{sf + K 4 (3) y( a ). 
AZf{s) = ( 4 ) x(,) 4 + K 4 (l) y( S ) 3 + i^ 4 (2) x( S )y( S ) + ^ 4 (3) y(s). 

4Zf(,)= ( 4 ) x( S ) 4 + ^K 4 (l) x(%( S K,( S ) + |^(2) y( S )z 2i0 ( S )+K 4 (3) y( a ) =AZf{s). 

4Z?(8) = Q x(,) 4 + K 4 (2) x(s)y(s) + 2^ 4 (3) y(s). 
AZf{s) = ( 4 ) x(,) 4 + K 4 (2) y(sf + 2K 4 (3) y(s). 
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4Zf( S ) = ( 4 ) x( S ) 4 + 3^(2) x(s)y(s) = 4Z™(s). 
AZf{s) = Q x( S ) 4 + K 4 (2) (2x(s)y(s) + y( S ) 2 ). 
AZf (s) = Q x( S ) 4 + K 4 (2) (x(s)y(s) + 2y{sf). 

AZf{ S ) = Q x( S ) 4 + K 4 (2) (x(s)y(s) + V2y(s)z 2 , (s)). 
4Z»( S ) = Q x(s) 4 + K 4 (2) (y( S ) 2 + ^yWza.oW). 
4Zf (,) = ( 4 ) x(s) 4 + K 4 (2) (x(s)y(s) + f + ^i/jW). 

F ~ Z 4 : 

4Z«(s) = Q x(s) 4 + K 4 (2) + *r p (B) x( a )(z lil/4 (a) + z 1>3/4 (s)). 

AZf{s) = (p x(s) 4 + K 4 (2) y(sY + tr p (B) (z 1 , 1/4 (,) 2 + ^(s) 2 )- 
4Z*°( S ) = Q x(s) 4 + K 4 (2) x{s)y{s) + ir p (B) (z 1)1/4 (s) + z lj3/4 (s)). 
AZf{s) = ( 4 ) x( S ) 4 + K 4 (2) y{sf + v^2 tr p (B) z 2 , 1/2 (s). 
F ~ Z 6 : 

6Zf( S ) = ( 4 ) x(s) 4 + %r p {B 2 ) (z 1A/3 (s)z 3fi (s) + z^/sO^oO*)) + K 4 (3) y( S ) + 
ir p (B) (^i,i/ 6 (s) +^i,5/e(s))- 

F ~ L> 3 : 

6Zf (a) = ( 4 ) x(s) 4 + 3^(2) yW^o W + %r p {B) (* 1>1/s ( s ) + ^WKoOO. 
F ~ Z| : 

8Zf ( S ) = ( 4 ) x( S ) 4 + 3^ 4 (1) x( S ) 2 y( S ) + F; 4 (2) (2x(s)y(s) + y(,) 2 ) + i^ 4 (3) y(s). 
8Z™(s) = ( 4 ) x( S ) 4 +K 4 (l) (2x( S ) 2 y( S ) + x( S )y( S ) 2 )+K 4 (2) (x(s)y(s) + 2y( S f)+K*(3) y(s). 
8Zf (s) = ( 4 ) x( S ) 4 + K${\) (x(s) 2 y(s) + x(s)y(s) 2 + y(sf) + 3K 4 (2) x( S )y( S ) + K 4 (3) y( S ). 
8Zf ( S ) = ( 4 ) x( S ) 4 + K 4 (l) (2x( S ) 2 y( S ) + x( S )y( S ) 2 ) + 3K 4 (2) x(s)y(s) + K 4 (3) y( S ). 
8Zf (5) = ( 4 ) x( S ) 4 + K*(l) (2x(,) 2 y( S ) + y( S ) 3 ) + K 4 (2) (2s(s)y( S ) + y(s) 2 ) + K 4 (3) 
8Z»(a) = 8Z p 34 ( S ). 

8Zf (5) = ( 4 ) x( S ) 4 + K 4 (l) (x( S ) 2 y( S ) + 2x( S )y( S ) 2 )+K 4 (2) (2x( S )y( S ) + y(,) 2 )+^ 4 (3) y( S ). 

8Z 40 ( S ) = 8Z» («). 
8Zf(s) = 8Zf 

8Z 42 ( S ) = ( 4 ) x( S ) 4 + 3^ 4 (1) x( S )y( S ) 2 + 3K 4 (2) *( s )y( s ) + i^ 4 (3) y( S ). 

8Zf (5) = ( 4 ) x( S ) 4 + K 4 (l) x(sfy(s) + K±{2) (2x(s)y(s) + y(s) 2 ) + 3^ 4 (3) y(s). 
8Z*%s) = ( 4 ) x( S ) 4 + K 4 (l) x(s)y(s) 2 + 3^ 4 (2) x( S )y( S ) + 3K 4 (3) y(s). 
F ~ Z 2 x Z 4 : 

8Z p 57 ( S ) = ( 4 ) x( S ) 4 + K 4 (l) x( S )y(s) 2 + K 4 (2) x( S )y( S ) + ir p (Bi) (z 1>1/4 (s) + z 1>3/4 (s)) + 
tr p (BiB 2 ) x( S )(z lil/4 (s) + z 1)3/4 (a)) + ^(3) y(a) = 8Z p 58 ( S ). 

F ~ F> 4 : 

8Z p 54 ( S ) = ( 4 ) x( S ) 4 + V2 F: 4 (l) x( S )y( S )z 2i0 ( S ) + F: 4 (2) (y( S ) 2 + y(s)z 2 , (s)) + 
y/2tr p (Bi) z 2jl/2 (s). 

8Z P 6 ( S ) = x ( s ) 4 + K p( 2 ) y(s)(x(s) + V2z 2fi (s)) + tr^) (* 1>1/4 ( S ) + z lj3/4 ( S )) + 
2F: 4 (3) 

8Z p 60 ( S ) = ( 4 ) x( S ) 4 + F; 4 (2) (3x(s)y(s) + V2y(s)z 2fi (s))+tr p (B 1 ) x(s)(z ljl/4 (s) + z li3 /4(«)). 
8Z p 61 (.) = ( 4 ) x(,) 4 + K 4 (2) (x(s)y( S ) + 2y(,) 2 + ^2 y(s)z 2 , (s)) + 

*r p (Bi) x(s)(«i,i/ 4 (s) + «i,s/4(a))- 
8Z p 62 (s) = ( 4 ) x( S ) 4 + K 4 (2) (3x( S )y( S ) + y(s)z 2 ,o(s)) +tr p (Bi) y(a)(«i,i/ 4 (s) + 2i,3/4(«))- 
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Remark 3.8. (a) With a similar asymptotic argument as the one used in Lemma 3.1, one can 
show that the functions x(s),y(s),Zi ! j ( (s) for 1 < i < 3, < k < \F\ — 1, appearing in the 
previous theorem are algebraically independent. This fact will be used in the proof of the main 
results in the following section. 

(b) Let r be a 4-dimensional Bieberbach group with canonical lattice and holonomy group F. 
We note that one can read off some properties of T from the expression of the zeta function. For 
instance: (i) V has diagonal holonomy representation if and only if Z p (s) 6 jpjZ[x(s), y(s)]. (ii) 

My is orientable if and only if for every i,j, k with i+j + k odd, the coefficient of x(s) l y(sy z(s) k 
in Zq (s) is zero. 

4. ISOSPECTRAL CLASSIFICATIONS 

In this section we shall use the information on the p-heat trace polynomials in the previous 
sections to determine all isospectral classes for the different type of spectra for all compact flat 
manifolds of dimension 4 in the class considered. We will thus find all pairs that are isospectral, 
p-isospectral, L-isospectral or [L] -isospectral in this case. 

Theorem 4.1 (p-isospectrality) . Let T,T' be two 4-dimensional Bieberbach groups with canon- 
ical translation lattice. Then: 

(i) If p = 0,4: Mp and Mr' are p-isospectral if and only i/{r,r'} is one of the following: 
{3,3"}, {3', 3"'}, {6,6'} (F ~ Z 2 ), {47,47'} (F ~ Z 3 ), {7', 10}, {8, 9'}, {14,14'}, 
{15,15'}, {18,19}, {19', 21}, {22,22'}, {24,26} (F ~ Z\), or {34,38}, {35,40}, 
{39, 41} (F ~ Zl), {57, 58} (F ~ Z 2 x Z 4 ). 

(ii) If p = 1,3: Mp and M-p> are p-isospectral if and only ifT,T' both belong to one of the 
following families: {3,3"}, {3', 3"'}, {5,5', 6, 6'} (F ~ Z 2 ), {47,47'} (F~Z 3 ), {7,9}, 
{7', 8, 9', 10}, {8', 10'}, {10", 11}, {14,14'}, {15,15'}, {18,19,20}, {19', 20', 21}, 
{22,22'}, {23,23'} (F ~ Z\), {24, 25, 26, 27, 28, 29, 29', 50, 51} (F ~ Z|,Z 4 ;, {34, 
36, 38}, {35, 40}, {39, 41} (F ~ Z\), {57, 58} (F ~ Z 2 x Z 4 ; or {60, 61} (F ~ D 4> J. 

fmj If p = 2: Mr and Mr' ore 2-isospectral if and only ifT,T' both belong to one of the 
following families: {2, 2', 2", 3, 3', 3", 3"', 4}, {6,6'} (F ~ Z 2 ), {47,47'} (F ~ Z 3 j, 
{7, 8, 9', 10', 10", 11', 12', 18, 19, 19', 21, 21', 23, 50}, {7', 8', 9, 10, 11, 12, 20, 20', 23', 
51} (F ~ Z^,Z 4 ;, {13, 13', 14, 14', 15, 15', 22, 22'}, {24,26} (F ~ Z^, {33,37,39, 
41, 43}, {34, 38}, {35, 36, 40, 42, 44} ~ Z|;, {57, 58} (F ~ Z 2 x Z 4y ). 

Proof. By Theorem 1.3, if Mr and Mr' are p-isospectral then \F\ = \F'\. Therefore when looking 
for p-isospectrality we only need to compare pairs both with F ~ Z 2 ,F ~ Z3, F ~ Z2 or Z 4 , 
F ~ Z\,lh 2 x Z 4 or D 4 or else F ~ Z@ or D3. Also, for each p, /3 P gives the multiplicity of the 
eigenvalue of A p , hence we need only consider pairs having the same (3 P . 

Take \F\ = 2. Here the groups are 2, 2', 2", 3, 3', 3", 3'", 4, 5, 5', 6,6'. Looking at the cor- 
responding p-heat trace polynomials in Section 3, we see that they are all different from each 
other, hence they can not be pairwise isospectral. However, the vanishing of the Krawtchouk 
polynomials for some values of p yields some p-isospectrality for 1 < p < 3. For instance, for 
p = 1 or 3 we know that K*(2) = 0, so Z*(s) = zf{s) = Z*(a) = Z*' (s) = 2x(s) 4 . Thus 5, 5', 
6 and 6' are l-isospectral and 3-isospectral. In the same way, Kp(l) = Kp(3) = for p = 2 and 
so Z 2 (s) = 3x(s) 4 for all groups 2, 2', 2", 3, 3', 3", 3"', 4, hence these groups are 2-isospectral to 
each other. 

In the cases of the remaining groups, i.e. those with \F\ = 3,4,6 and 8, if p = 0, using 
the complete list of the p-heat trace polynomials given in Section 3, we easily see that only 
the pairs {3,3"}, {3', 3"'}, {6,6'}, {47,47'}, {7', 10}, {8,9'}, {14,14'}, {15,15'}, {18,19}, 
{19', 21}, {22,22'}, {24,26}, {34,38}, {35,40}, {39,41}, {57,58} have the same p-heat 
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trace expressions for p = 0. Actually they have the same p-heat trace for all p, thus they are 
p-isospectral for < p < 4. 

In the case of the p-spectrum for 1 < p < 3, again by comparison of the p-heat traces one 
obtains the p-isospectral sets asserted in the theorem. □ 

Corollary 4.2. Let T,T' be as in Theorem 4-1- Then Mr and Mr' are isospectral if and only 
if they are p-isospectral for < p < 4. 

Remark 4.3. (i) Observe that the groups in each isospectral pair in Theorem 4.1 have the same 
holonomy representation. 

(ii) The assertion in the corollary fails to be true for general flat manifolds. Examples of 
0-isospectral and not 1-isospectral flat manifolds in dimension n > 6 are given in [6]. In [2] two 
nonhomeomorphic flat 3-manifolds that are isospectral on functions are given, showing this is 
the only such pair. These manifolds have different holonomy groups: Z4 and Z2 x 7h<i respectively 
and one can check they are not isospectral on 1-forms by comparison of the p-heat traces. 

(iii) The theorem also shows that isospectrality on functions is much less common than p- 
isospectrality for 1 < p < 3. Some facts can be observed. For p = 1,3, one of the p-isospectral 
pairs, {60,61}, involves two groups with non-abelian holonomy group D4. Also, in the p- 
isospectral set {24, 25, 26, 27, 28, 29, 29', 50, 51} five different holonomy representations are 
present. Here, the first seven groups have holonomy F ~ 1^ and give orientable manifolds, while 
the last two groups have holonomy F ~ Z4 and produce non-orientable manifolds. Relative to 
2-isospectral classes, again some manifolds are orientable and others not, and different holonomy 
groups occur. 

We now turn into the study of length spectra for the flat 4-manifolds considered in this paper. 
It is a classical result that in the case of flat tori the eigenvalue spectrum and the length spectrum 
determine each other. This is not true for general Riemannian manifolds, even in the flat case. 
There exist 4-dimensional flat manifolds (namely 8,9' in our notation) which are isospectral but 
such that some multiplicities of closed geodesies are different (see [8], Ex. 3.4, up to isometry; 
see also the references therein for earlier examples in the case of nilmanifolds) . In particular, 
this shows that in general, the eigenvalue spectrum does not determine the length spectrum. In 
the converse direction, in [8] an example was given, for n = 13, of [L] -isospectral flat manifolds 
which are not isospectral on functions. 

The study of [L] -isospectrality is delicate since multiplicities must be taken into account. We 
recall that the multiplicity of a length / is the number of conjugacy classes in T having length 

We shall see that already in the context in this paper, most L-isospectral pairs fail to be 
[L]-isospectral. The next result describes the situation. 

Theorem 4.4 (L-isospectrality). Let T,T' be two 4-dimensional Bieberbach groups with canon- 
ical translation lattice. 

(i) The L-isospectral sets are the following: {3,3"}, {3', 3"'}, {6,6'}, {47,47'}, {7', 10}, 
{8, 9'}, {14, 14'}, {15, 15'}, {18, 19}, {19', 21}, {22, 22'}, {24, 26}, {34, 38, 39, 41}, {35, 40}, 
{57,58} and, in addition, the sets {25,27}, {33,43}, {42,44}. 

(ii) The \L\-isospectral pairs are {3,3"}, {3', 3'"}, {6,6'}, {47,47'}, {7', 10}, {14,14'}, 
{15,15'}, {18,19}, {19', 21}, {22,22'}, {24,26}, and {57,58}. 

Proof, (i) All isospectral pairs in (i) of Theorem 4.1 are L-isospectral. As we have seen in 
Section 3, the information on the L-spectrum of T depends only on the monomials ocurring 
in the 0-heat trace. Since for the L-spectrum we do not look at multiplicities, T,r' will be 
L-isospectral if and only if the same monomials occur in the 0-heat trace polynomials (see 
corollary 3.6). 
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Checking the list of p-he&t traces we see, for instance, that 4Zp 5 (s) = ( ) x(s) 4 +Kp (2) (y(s) 2 + 
2x{s)y(s)) and4Z2 7 ( s ) = (J i(s) 4 +K p 4 (2) (x(s)y(s)+2 y(s) 2 ) , thus 25 and 27 are L-isospectral. 
In the same way we see that {33,43}, {42,44} and {34,38,39,41} are L-isospectral. 

(ii) Since [L]-isospectral implies L-isospectral, we must only consider the sets given in (i). 

In [8], the length spectrum of flat manifolds is studied and a criterion to prove [L]-isospectrality 
is given (see Proposition 3.1). This criterion depends on the existence of a bijection with cer- 
tain properties between partitions V and V' of F and F' respectively. Also, in [8] the criterion 
is applied to show [L]-isospectrality of some pairs, including the pair {24, 26} (see Example 3.3, 
up to isomorphism). 

We shall sketch the verification for the pair {3, 3"} (see Table 2.1). We have that V (resp. V") is 
generated by 7 = d(I, J)L b (resp. 7" = d(I, J)L b ») and L A , where b = ^, (resp. b" = 6l+e 2 3+e4 ). 

The elements of V (resp. T") are either of the form 7L A (resp. j"L\) or L\ with A = m\ei + 
m2e2 + m^es + 772464 G A, i.e. rrii G %. 

We define a correspondence <f> from T to T" such that L\ ^ L\ and ^ 7"L A L_ e3 . We 
note that both of these elements have squared length (mi + i) +7712 + ^("1-3 + TO4) 2 . 

If we use the relation (4.1) below we see that 7L A ~ 7-^A-^m(e 3 -e 4 ) an d similarly 7"L A ~ 
7" L\L m ( e3 _ e4 -j, for any m G Z. We note that (4.2) does not introduce new conjugacy relations 
among the elements of T and T" , since the holonomy group is cyclic. 

It is now easy to see that the map 4> induces the bijection & = Id from F to F" = F, satisfying 
all the conditions in Proposition 3.1 of [8], hence {3,3"} are [L]-isospectral. 

The [L]-isospectrality of the remaining pairs in (ii) of the theorem, can be proved similarly, 
with = Id as the bijection. We shall omit this verification. 

To conclude the proof of the theorem we must show that the remaining L-isospectral sets, 
{8,9'}, {25,27}, {33,43}, {34,38,39,41}, {35,40}, {42,44}, are not [L]-isospectral. 

The pair {8, 9'} was shown not to be [L]-isospectral in [8] (up to isometry, see Example 3.4). 
We now discuss in some detail the case {25, 27}. 

If / is the length of a closed geodesic in Mr, the multiplicity is given by mr(I) = 2^ 7 eA\r m 7(0 
where 

m 7 (0 = #{[ 7 L A ] G [r] : A G A, £( 7 L A ) = /}. 
In order to compute these multiplicities it is necessary to have a parametrization of the T- 
conjugacy classes in T. This is complicated in general but it simplifies for Bieberbach groups 
of diagonal type. From [8] we have the following relations. For general T, if 73 = BiL^^j = 
BjLfr. G T and A, fi G A,the conjugations 7iL A 7j _1 and L M (7jL A )L_ M give the relations: 

(4-1) L\ ~ L BiX and 7i L A ~ liL\L {B .-i_ Id)ll 

Furthermore, if F is abelian, the conjugation jj^Lx^j -1 , with i / j, gives the relation 

(4.2) 7i L A ~ JiL VjAX) where rj^(A) := BjX + (Bj - Id)h + B, (B { - Id)bj 

The next table gives the elements 71, = L?i 61,72 = L ? 2L&2,7i2 = Bvibvi in 25, 27, in column 
notation, i.e. showing in the columns, for each BL b G T, the diagonal entries of B together with 
the coordinates of the corresponding translation vector b (see [6]). 



Bx bf bf 


B 2 bf bf 


-BiL? 2 6^2 ^12 


-1 
-1 

1 1/2 
1 1/2 


1 1/2 1/2 
-1 1/2 1/2 
-1 

1 1/2 


-1 1/2 1/2 
1 1/2 1/2 

-1 1/2 
1 1/2 1/2 



The squared lengths of closed geodesies corresponding to elements not in A are given by 
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25: m§ + (m 4 + ±) 2 , (mi + ±) 2 + m|, (m 2 + J) 2 + (m 4 + \f , 
27: (m 3 + ±) 2 + m 2 , (mi + ±) 2 + (m 4 + \) 2 , (m 2 + |) 2 + (m 4 + |) 2 , 

with mj 6 Z. 

The minimal length in both cases is / = ^. We shall show that this length has different 
multiplicies for 25 and 27. In 25 the elements with I = \ have the form: 

B\L24 +Xl where Ai = m\e\ + m 2 e 2 — m^e^ with mi,m 2 G Z, m 4 G {0, 1} 

B 2 Lei +X2 where A 2 = — niei + n 2 e 2 + n 3 e 3 with rc 2 , n 3 G Z, ni G {0, 1}. 

In 27 the elements with / = \ have the form: 

B\L^_ +Xi i where X\ = nei + r 2 e 2 — r 3 e 3 with n, r 2 G Z, r 3 G {0, 1}. 

It follows from (4.1) that mi, m 2 , n 2 , n 3 , n, r 2 can be taken mod (2). 
Using the relations (4.2) we have: 

7l^miei+m 2 e 2 +m4e4 ~ 7l-^(mi-l)ei-(m 2 +l)e 2 +m4e4 m ^5 

72-^niei+n 2 e 2 +n 3 e 3 ~ 72-^-(ni+l)ei-(n 2 +l)e 2 +n 3 e3 m 25 

7l^riei+r 2 e 2 +r 3 e3 ~ 7l-^(ri-l)ei-(r 2 +l)e 2 -r 3 e 3 m 27. 

These three relations divide by 2 the number of relations, and there are no other relations. 
Thus m25(^) =4 + 4 = 8 while 777-27(7?) = 4. This shows that 25,27 do not have the same 
multiplicities. 

The verification of non [L]-isospectrality in the remaining cases can be done similarly as for 
25, 27 (by comparing the multiplicities of small lengths) and will be omitted. □ 
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